Bounds to unitary evolution 
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Upper and lower bounds are established for the survival probability \(tp(0)\ip(t)}\ 2 of a quantum 
state, in terms of the energy moments (tjj(0)\H n \ip(0)) . Introducing a cut-off in the energy generally 
enables considerable improvement in these bounds and allows the method to be used where the 
exact energy moments do not exist. 
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I. INTRODUCTION 

How rapidly (or how slowly) can a state evolve? This 
question has been approached [l|, 0, H, 01 through the 
survival probability of the state: 

P(t) := McM-*Ht/h)i;)\ 2 = |(exp(-zm/^)| 2 , (1) 

which is the probability that the system, initially in state 
ip, will be found to be still in that state after time t. The 
Hamiltonian H is assumed to be independent of the time. 
Lower bounds to P(t) have been established [l|, y] in 
terms of the energy uncertainty AE := {{H — (iJ)) 2 ) 1 / 2 , 
but these bounds can be well below the actual evolution. 
Also it has been claimed @ that there can be no upper 
bound (other than P(t) ^ 1) in terms of AE alone; with- 
out an upper bound one can not be sure that the state 
will change at all. Here, upper and lower bounds will be 
found in terms of higher energy moments. 

Expanding the exponential in Eq. ([T]) as a power series 
gives 



P(t) = l~(h 2 -h 2 1 ) — + (h i -4h 3 h 1 +3h 2 2 ) 
nr 



I2h 4 



(2) 



where h n := (H n ), the n-th energy moment. The absence 
of a linear term in this expansion is important to the 
discussion of the quantum Zeno effect d, 0] • 

There are good reasons why hi and hi should exist [1] , 
but for many states used in physics some of the higher 
energy moments do not exist. It will be established in 
Section HTT1 that, apart from the factor (—1)™, the coeffi- 
cient of t 2n in Eq.([2]) must be positive if it exists, and that 
the partial sums of this series give alternately upper and 
lower bounds to P(t); an upper bound if the last included 
term is positive or a lower bound if the last included term 
is negative. Section ITVl introduces a cut-off in the energy, 
equivalent to projecting onto a finite-energy subspace. 
This enables considerable improvement in these bounds 
(including the one in terms of AE) and also allows the 
method to be used where the exact energy moments do 
not exist. We first need to show that the coefficients in 
Eq.® are moments over the autocorrelation of the en- 
ergy distribution. 
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II. THE ENERGY DISTRIBUTION 

Consider a complete set of commuting observables 
H, K with common eigenstates \4>e,k), so that H\(f>E, K ) = 
E\<Pe,k) and K\<pE, K ) = k\<Pe,k)- (In general, K repre- 
sents a set of operators and k a set of eigenvalues.) Then, 
for any function f{H) of the Hamiltonian, 



(f(H)) = dE dKf(E)\(<f> E , K 



(3) 



Define the energy distribution p(E) 
Then J dE p{E) = 1 and 



/^|(^,k|^)| 2 



(f(H)) = I dEp(E)f(E). 



(4) 



For example, h n := (H n ) = J dE p(E) E n . 

Let L and M be the lower and upper bounds to the 
energies for which p{E) is non-zero, but allow the pos- 
sibility that M = oo and even the unphysical case that 
L = — oo. The survival probability P(t) is not changed 
by a shift in energy and for all the examples used here 
L = or — oo. In terms of p(E), 



P(t) 



? i\ME' -E)t/h 



M nM 

dE / dE'p(E)p{E')e l 

L J L 

dE / dE'p{E)p{E') cos 

L JL 



E 



4(5) 



Change the integration variables from {E, E'} to {E, e} 
with e := E' — E, and introduce the autocorrelation of 
p{E) through the even function 



W(e) 



M-e 



dE p(E)p(E + e) fore > 0. (6) 



Then W(e) is never negative and 



Pit) = I 
Jo 



M—L 



deW{e) cos( 



e t , 



(7) 



Also / M L deW{e) = P(0) = 1. 

Expanding the cosine in Eq.©, gives 



P(t) 



e 2 t z 
2\h 2 



e 4 t 4 
4!F 



e 6 t 6 
6lh 6 



(8) 



2 



[which must agree with Eq. (J2J)] , where 

•M-L 

deW{e)e r> 



(9) 



The moments e n are positive and can be expressed, for 
even n, in terms of the energy moments hk with k < n : 



M-L 



L—M 



deW{e)e n = J J de dE p(E)p(E + e)e" 
= JJdEdE' p(E)p(E')(E-E') n . (10) 

Now expand (E — E') n as a sum of products of powers of 
E and E'. Thus, for example, (25 - E') 2 = E 2 - 2EE' + 
E' 2 leads to e 2 = h 2 - 2h\ + h 2 = 2{h 2 - h\) = 2{AE) 2 
and e4 = 2(/i4 — 4/13/11) + 6/i| . More generally, 



n/2-l 

^ = 2 £ (- 

fc=0 



h n - 



k\ {n-k)\ 



iW2/ 



+ (-l) n/ ^( 



(n/2)! 



(11) 



1.0 




0.7 



FIG. 1: Bounds to the survival probability for the energy 
distribution p(E) oc E~ 1/2 e~ B/l , E > 0. The time unit is 
ft/7. The dashed curve is the exact survival probability P(t) 
and the solid curves are the bounds given by Eqs. (|12ll4[l . 
The numbers on these bounds give the largest energy moment 
used. 



III. BOUNDS ON THE SURVIVAL 

We first show that the partial sums of the Taylor series 
for cos a; provide alternately upper and lower bounds to 
cos x. To prove this, note that if we have an upper (or 
lower) bound to cos a; for all x > 0, applying sinx = 
J cosudu gives an upper (lower) bound to sinx. Then 
applying cos x = 1 — JJf sin u du gives a lower (upper) 
bound to cosx. Thus cosx < 1 =>■ sin a; < x =>■ cos a; > 
1 — \x 2 sinx > x — -^x 3 cosx < 1 — + ^i 4 
and so on. Using this sequence of inequalities for cosx, 
it follows directly from Eq.([7} that 



1 . 




0.7 



P(t) 
P(t) 
P{t) 



> 1 



< 1 



> 1 



e 2 t z 
2h 2 
e 2 t 2 
2\h 2 
e 2 t 2 
2\h 2 



= 1 - (AE) : 



e 4 t 4 
4!fi 4 

4\h 4 



e 6 t 6 
6lh 6 



(12) 
(13) 
(14) 



and so on. Whereas the series Eq.© for P(t) may or 
may not converge, each of these bounds is valid provided 
the moments e„ in it exist. 

The lower bound in Eq. (|12p is well known and has been 
improved [l], 9 to P(t) > cos 2 (A£ t/K). To my knowl- 
edge, the other bounds are new; furthermore I know of 
no other upper bounds to P(t). 

As a simple example, consider the energy distribu- 
tion p(E) (x E~ 1 / 2 e- E /i, E > 0. [One of many pos- 
sible realizations of this as a wavefunction is as the free 
Gaussian if) = (t — it)^ 1 / 2 exp( j^imx 2 /(t — it)), where 
r = 5 ft/7-] The exact evolution has the survival prob- 
ability P(t) = (1 + {jt/h) 2 )- 1 / 2 and the coefficient of 
t n in the power series for this must be e n /n\. In fact, 
e n = (n — l) 2 (7/fi,) 2 e„_2 with eo = 1. The autocorre- 
lation is W(e) = 2 / (ttj) K (e / j) , but it is not needed. 



FIG. 2: Bounds to the survival probability for the energy dis- 
tribution p(E) oc (1 + E/-yy 7/2 , E > 0. The upper curve 
is the exact survival probability P(t) and the solid curve la- 
belled "2" shows the bound given by Eq, (|12|l , The dashed 
curve is cos 2 (AEt/h), which is the best possible lower bound 
if only AE is known. The time unit is ft/7. 



All energy moments exist because the energy distribu- 
tion falls off faster than any power of the energy. Fig. [1] 
shows the exact evolution and four successive bounds for 
this case. 

When the energy distribution p(E) decreases slowly 
(slower than exponentially) then higher energy moments 
may not exist and even if some do exist the bounds pro- 
vided by Eqs. (|12j|14")) may be very poor. For example, if 
p{E) oc (1 + E/-1)- 7 / 2 , E > 0, then e 2 = f 7 2 but no 
other even moments exist. Eqs. (|12fll4[) provide no upper 
bound and the lower bound P(t) < 1 — \e 2 t 2 /h 2 is very 
poor, as shown in Fig. [5] We will now show how using 
an energy cut-off yields good upper and lower bounds for 
this system and many others. 
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IV. USING A CUT-OFF IN THE ENERGY 

Following the work of Uffink and Hilgevoord @, we 
cut off the energy at say E — c and write 



(15) 



1 . 



p{E) dE. 

The exact state of the system can be expressed as 

^(t) = V5$(t) + vT=^x(*). 

where, in the notation of Eq.Q, 



-LjdnJ dEe- lEt / h p(E)<t> E , K , 



(16) 



(17) 



Xit) := 



A I 



dEe- lEt ' h p{E)4> E ^. (18) 



Then ip(t) and x(t) are normalized and orthogonal, and 
the survival amplitude A(t) := (ip(0)\ip(i)) is 



A(t) = aA(t) + (1 - a)B(t), 



(19) 



where Ait) := (ip(0)\ip(t)) and := (x(0)|x(i))- Since 
\B(t)\ < 1 we have the inequalities [8j 



\A(t)\ < a\A(t)\ + (1 - a) 
\A(t)\ > a\A{t)\-(l-a). 

We now apply the bounds in Eqs.(10)-(12) to \A\; 

e 2 t 2 



(20) 
(21) 



\A(t)\ > a (l-^/ 2 -(l-a) 



\A(t)\ < a(l- 

\A(t)\ > a(l- 
and so on, where 



e 2 t 2 

m 2 

21H 2 



e 4 t 4 
4!F 



) 1 / 2 + (l 
e e t 6 



a) 
) 1/2 "(1 



c-L 



deW(e)e n , 

e 

dEp(E)p{E + e). 



(22) 
(23) 
a) (24) 

(25) 
(26) 



Note that it is not necessary to calculate the autocorre- 
lation W for the truncated system; the moments e„ can 
be obtained from the energy moments of the truncated 
system using Eq. fTTj) . 

Returning to the example with p(E) oc (1 + E/^)~ 7 / 2 
which has a = 1 — (1 + c/j)~ 5 / 2 , Fig. [3] shows these 
bounds for selected values of c using the quadratic lower 
bound in Eq.(23) and the quartic upper bound in Eq.(24). 

It is now clear that the bounds from a given cut-off c 
are good only for a limited range of the time and the best 
we can do with each of the bounds in Eqs. (|22ll24"|) is to 
calculate the envelope as c varies. Each bound to \A(t)\ 
is given by a function 



y(t, c) := a{c)\Jp n {t,c) ± (1 - a(c)), 




0.4 



FIG. 3: Bounds to the survival magnitude for the energy 
distribution p{E) oc (1 + E/^)~ 7/2 , E > 0. The dashed 
curve is the exact magnitude |A(i)| of the survival and the 
solid curves are examples of the bounds in Eqs. (|22j | !£?fl for 
selected values of the cut-off energy. The lower bounds come 
from the quadratic form in Eq, (|22[l and the upper from the 
quartic form in Eq. (|23p . The numbers on these curves give 
the values of c in units of 7. The curve with c = 00, i.e. 
no cut-off, corresponds to the quadratic form in Eq. (|12p ; this 
also appears in Fig. [2] 
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\A(t)\ 
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FIG. 4: Bounds for the same system as in Fig. [3] The dotted 
curves are the envelopes of the bounds given by Eqs. (|22l24[l , 
The numbers on these envelopes give the largest energy mo- 
ment used. 



where p n {t, c) is a polynomial of degree n in t. The en- 
velope of these bounds as c varies is found by solving 



d c y(t,c) = 



(28) 



(27) 



to give t as a function of c. This can be done explicitly 
for the quadratic case, n = 2, and for the quartic, n = 4; 
but in general only numerical solution is practical. The 
details are in the appendices. Then inserting this i(c) 
into y(t, c) gives the envelope y(t(c), c) at time t{c) para- 
metrically in terms of c. The envelopes for the first four 
bounds for the distribution p(E) oc (1 + E/~f)~ 7/2 , E > 
are shown in Fig. [3J 
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FIG. 5: Bounds for the Breit-Wigner distribution. The solid 

of the survival, 



curve is the exact magnitude |A(t)| = e -7 *^ 
and the dotted curves are the envelopes of the bounds given 
by Eas. (|22j | 24]) . The numbers on these bounding curves give 
the largest energy moment used. The time units are ft/7. 
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FIG. 6: Bounds for the real part and the negative of the imag- 
inary part of the survival amplitude A(t) for the distribution 
p(E) oc E~ 1/2 e~ Eh , E > 0. The dashed curves show the 
exact amplitude and the solid curves are the bounds given by 
Eqs. (|29l32p . The numbers on these curves give the largest 
energy moment used. 



The Breit-Wigner system p(E) oc [(E - E ) 2 + 7 2 ]- 1 
is unphysical because neither (H) nor (H 2 ) exist, but 
it is used because it exhibits exact exponential decay: 



\A(t)\ 



-~/t/h 



None of the bounds in Eqs. (|12IU"4")) can 



be used, but cutting off all energies outside the range 
Eq — c to Eq + c gives the bounds shown in Fig. [5j 



V. BOUNDS ON THE REAL AND IMAGINARY 
PARTS OF THE SURVIVAL AMPLITUDE 



The survival amplitude can be expressed as A(t) = 
f dEp(E)exp(-iEt/h) = R(t) - il\t). The methods 



used in Section [HT] applied to R(t) and I(t) lead to 
hit 



I(t) 
R(t) 
I(t) 



< 



> 1 



M 2 

2ft 2 



> 



hit 



R(t) < 1- 



M 3 
3!ft 3 
h 2 t 2 hA 4 



2!ft 2 4!ft 4: 



(29) 
(30) 
(31) 
(32) 



and so on. Fig. [5] shows the results when this is applied 
to the same distribution as used in Fig. [TJ Again an 
energy cut-off could be used to improve these bounds, or 
to apply them when the energy moments do not exist. 



VI. DISCUSSION 

All cases considered here have just one t(c) > from 
Eq. (|A2p for each n, but this has not been proven in gen- 
eral. Note that if p{E) increases suddenly then t(c) may 
not be monotonic in c; then there is more than one bound 
for a period of time and the only the best one is of value. 

Thanks are due to M J W Hall for useful comments. 



APPENDIX A: DETERMINING THE ENVELOPE 

Applying d c y(t,c) = to y{t,c) in Eq. ([27]) leads to 
Pn(t,c) = ip n (t,c) + £rd cPn {t,c)) 2 . From Eq.©, 
d c e n = — 2(a'/a)(e„ — b n ), where a' := d c a = p{c), 



b n :=- fdEp(E)(c-E) n = V 



{-l) k n\ 
k\{n-k)\ 



(Al) 

and ftfe := a 1 j L dE p(E)E n . Then the equation for the 
envelope is 



EM) 



fc/2 



e k t k 

h k k\ 



( E 



b k t k 



h k kl 



(A2) 



This must be solved for tic) > and the envelope is 
p n (t(c),c) at time t(c). 



Replacing E-E' by (c-E')-(c-E) in Eq.^ shows 
that e n /nl = J2k=o(~ 1 ) k B kB n -k, where B k := b k /kl; so 
the envelope equation is simplest in terms of the only. 
For the quadratic bound (n = 2): t(c) = 2/161/62. The 
quartic bound (n = 4) leads to a cubic equation in t 2 : 
B 2 -2BiB 3 f]+2B 2 B 4 r] 2 -Bj7 1 3 = with 77 := t 2 /h 2 . The 
relevant solution of this is t(c) 2 



where d\ := 6163 



bt do 



8h 2 (b 2 -d 2 + di/d 2 )/b i , 
d 3 ) 1/3 , with 



d 2 fl 2 



2, 

(16&H - 246i6 2 6 3 + 96 2 6 4 )/16. Note that all the 



b k > 0, and d\ > from the Schwarz inequality. Also 



16d 3 = (4&2 /2 - 36i64^) 2 + 246i6 2 (6 2 / "64 / " - 63) and this 
is positive because 6264 > 6 2 from the Schwarz inequality. 
It follows that this root of the cubic is positive. The 
next bound (n = 6) leads to a quintic equation in t 2 and 
numerical solution is probably the most practical option. 



(1/2 



1/2,1/2 
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FIG. 7: Bounds for the survival magnitude |A(t)| for p(E) — 
1/M, < E ^ M. The solid curve show the exact amplitude, 
the dashed curves are the bounds for n = 2, 4, 6, 8 and the 
dotted curves are the corresponding envelopes for varying cut- 
offs. The dotted curve labeled oo is the limit of the upper 
envelopes asn-> oo. The unit of time is h/M. 



APPENDIX B: FINITE RANGE OF ENERGIES 

When p(E) = for all E > M, the energy cut-off must 
also stop at c = M. This causes the envelope to be valid 
only for times greater than t(M) given by d c y = at 
c = M. But the envelope will match smoothly to the 
bound without cut-off at that time, because the enve- 
lope osculates the sequence of bounds as c approaches M. 
These two bounds together provide a continuous bound 
for all times until the envelope reaches either or 1. 

To illustrate this, consider the simple case of a square 
energy distribution: p(E) = 1/M for s$ E ^ M. With 
an energy cut-off at E = c, e n /n\ = 2c n /(n + 2)! and 
b n /nl — c n /(n + 1)!, so the envelope equation (|A2[) be- 
comes independent of c if t(c) — fir„/c, where r„ is a 
dimensionless constant. Each upper envelope has the 
form 1 — (1 — o n )t n /t and each lower envelope the form 
(1 + cr n )t n /t — 1, where r„ is the positive solution of 

" (-l) fc / 2 2r* » (-l) fc / 2 r* ,2 

^ (fc + 2)! - ( ^ (jfc + 1)! > ' (B1) 

even k even k 

t n := Hr n /M and c 2 equals either side of Eq. ljBip . 

The envelope for n = 2 has r 2 = 3 and a 2 = 1/2, 
which gives the bound to \A(t)\ as y 2 = 9H/(2Mt) — 1 
for t > 3H/M. Since y 2 = at Mt/H = 9/2, the useful 



range of this envelope is 3 < Mt/h < 4i. This can 
be seen in Fig. \7\ which also shows the bounds for n = 
4,6,8. The limit as n — > oo can be taken for the upper 
bounds by summing the series in Eq. ljBip in terms of 
sinr n and cosr„, leading to Too = 27r, Coo = and the 
bound \A(t)\ < 1 - 2-Kh/{Mt) for t > 2wh/M. 




FIG. 8: Bounds for the survival magnitude |j4(£)| for p(E) = 
0.7 5(E) + 0.2 5{E - M/2) + 0.1 5{E - M). The solid curve 
shows the exact amplitude, the dashed curves labeled 2 and 4 
are the bounds without cut-off, the dotted curves labeled 2c 
and 4c are the bounds with M/2 < c < M, and the line for 
c < M/2 is also shown. The unit of time is h/M. 



APPENDIX C: DISCRETE ENERGY SPECTRA 

The envelope does not exist in regions where the energy 
spectrum is discrete, i.e. where the energy distribution 
consists of (5-functions only, because the bounds provided 
by Eqs.([22" |) -(f2"4 |) do not change as c moves from one 8- 
f unction to the next. Then the present method gives 
a continuous series of bounds and t(c) can be used to 
specify the period of time that each bound will be valid. 

As a simple example, consider a 3-state system with 
equal space between the energy levels: p(E) — a S(E) + 
a x 8(E - M/2) + a 2 S{E - M) with a + a x + a 2 = 1. 
For < c < M/2, a — a , all the are zero, and the 
only bound that can be obtained from Eqs.([22|)-([24l is 
\A(t)\ ^ l-2o Q (useful only if a Q < |). If M/2 <c< M, 
a = ao + Oi, hk = (M/2) k ai/a and the bound from 
Eq.([22| is valid for i(M) < t < t(M/2), where t(c) = 
2hb 1 /b 2 , which leads to t(M/2) = Ah/M and t(M) = 
(H/M)(2a + oi)/(a + |oi). For t < t(M), one must 
use the bound without cut-off. This is illustrated for a 
particular choice of the ak in Fig.0 which also shows the 
quartic bounds. 
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